HANDS ON - LOW LEVEL 


Flights of fancy 


Any would-be Top Gun applicants roll up — this month, Mike Liardet presents a highly 
realistic and involving flight simulation program, guaranteed to sort the best from the rest. 


ite this month’s Low Level we will 
look at how the underlying model for 
a flight simulator can be built in Visual 
Basic. We will consider the simulation 
of a small single-engined trainer aero- 
plane such as the two-seater Cessna 
Aerobat (specification given in Fig 1). 
In many respects it’s more challenging 
to implement this type of plane in a 
simulation. High-performance jet fight- 
ers may be more difficult to fly in real- 
ity, but with such an excess of thrust 
compared with all the other forces act- 
ing on them, they can fairly accurately 
be modelled simply to fly wherever 
they’re pointing. 

The basic technique is outlined in 
the flow diagram in Fig 2. After ini- 
tialising the various variables and arrays 
(which may, for example, place the 
plane stationary on the ground, or in 
cruise at altitude) the program contin- 
uously executes the various stages in 
the main loop until either the plane 
crashes, or the user decides to aban- 
don the simulation. 

The main loop is executed in realtime, 


Fig 1 Cessna Aerobat: 
Technical Specification 
Engine: Rolls Royce O-240-A 130 HP 
(97 kilowatts) 

Fuel Capacity: 98 litres (85 litres 
usable) 

Empty weight: 520kg 

Maximum all-up weight: 750kg 
Top Speed at sea level: 113 knots 
Cruise at 75% power at 7000 feet: 
111 knots 

Range at cruise (no reserve): 

348 nautical miles 

Rate of climb at sea level: 

846 feet per minute 

Never Exceed speed: 168 knots 
Stall speed: 48.5 knots 


so in one second on the user’s watch the 
plane goes through the same changes in 
the simulation that a real plane would 
go through in one second. Of course, 
different computers have different 
speeds of execution, but in theory this 
doesn’t greatly affect the accuracy of 
the simulation — a fast computer sim- 
ply allows the main loop to be executed 
more times per second. 

Notice that one of the points in the 
loop updates the screen display, so with 
a fast computer processing the loop 
more quickly, the display gets updated 
more often. In contrast, a slow com- 
puter may be able to update the screen 
only a few times per second, leading to 
arather jerky presentation. This can be 
annoying for the user, but strictly speak- 
ing the computer speed doesn’t greatly 
affect the simulation in any other way. 
There can, however, be other problems 
with different processing speeds as the 
varying response times can affect the 
speed at which the simulator reacts to 
inputs from the user. This results in 
instability, causing the plane to go out 
of control. This is quite realistic, as 
somewhat alarmingly, real planes can 
suffer from this problem as well. 

The first step in the loop updates the 
screen display to reflect the plane’s cur- 
rent position. Most flight simulators 
show an instrument panel on the 
screen, possibly with a view of the out- 
side world given through a simulated 
windscreen. As a programming task it’s 
not too difficult to implement the 
instrument panel, simply by arranging 
to update the values given on the 
altimeter, air speed indicator and other 
dials. The instruments we might want 
to make available are shown in Fig 3. 

Unfortunately, providing a display 
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routine with anything that can provide 
areasonable picture of what’s going on 
outside the cockpit is a major under- 
taking. The simulated outside world 
first has to be modelled in some fash- 
ion, with the position, shape, texture 
and colour of objects, such as moun- 
tains, lakes, runways and buildings all 
set up in an appropriate data structure. 
There must be a display routine which 
can take this model in conjunction with 
the plane’s position and orientation in 
order to decide what the pilot would 
see through the windscreen, and then 
draw it. This involves some heavy pro- 
cessing, not to mention some fairly 
hairy mathematics, so we won’t take 
this any further here. 

Any adventurous programmers who 
would like to try should be warned that 
a Visual Basic implementation is 
unlikely to be fast enough for a real- 
time simulation. If it takes several sec- 
onds to draw the picture of the outside 
world, then obviously it will take even 
longer than this to go round the main 
simulation loop, and the simulation 
will be too unresponsive to be of any 
entertainment value. 

The next step in the loop looks at 
the user’s actions with the simulated 
controls. The user can provide some 
input here by changing the position of 
the plane’s joystick or control column 
and changing throttle settings. Moving 
the joystick to the left or right activates 
the ailerons. These are hinged flaps 
mounted on the trailing edge of the 
wings, and working in opposition to 
each other. Moving the joystick to the 
left causes the plane to roll to the left 
and likewise moving it to the right 
causes a right roll. When using the 
ailerons the plane keeps on rolling, 
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Yes 


Fig 2 Flow diagram for the flight simulator 


possibly even going inverted, until the 
ailerons are neutralised. 

Moving the joystick backwards and 
forwards activates the elevators. These 
are hinged flaps mounted on the tail- 
plane. Pushing the joystick forward 
pushes the elevators down, which tends 
to kick the tail-plane upwards causing 
the aeroplane to pitch downwards. 
Pulling the joystick back has the oppo- 
site effect, causing a pitch upwards. 
The elevators are generally used to 
balance the plane correctly for dif- 
ferent flying speeds. Unlike ailerons, 
a continuous elevator input need not 
necessarily cause the plane’s attitude 
to change through 360° if it is held 
indefinitely. 

Real aeroplanes also have a rudder 
which is mounted on the trailing edge 
of the fin and controlled by the pilot’s 
feet. Most people tend to imagine that 
it is used to turn the aeroplane in flight, 
but this is incorrect. It is used on the 
ground to steer the aeroplane and also 
in the air for aerobatic manoeuvres, but 


for ordinary flight it should only ever 
be used in conjunction with the 
ailerons to counteract a phenomenon 
called ‘adverse yaw’. An ‘ideal’ aero- 
plane, such as the one we are simulat- 
ing here, need not have any adverse 
yaw and so it does not need a rudder. 
This begs the question of how the 
plane’s heading is changed. 

The answer is that banking the plane 
(using ailerons) tilts the lift force on 
the wings so that in part the plane is 
pulled horizontally by the lift. This 
pull acts on the wings which are at the 
front of the plane, and the rear is effec- 
tively just dragged along behind, result- 
ing in a change of heading. In this con- 
text the fin and tail-plane act like the 
feathers on a dart, tending to point the 
plane in the direction of the airstream. 

The third step in the loop takes the 
plane’s velocity through the air and 
uses it to calculate all the forces acting 
on it. There are lift and drag forces from 
the wings, tail and fin, plus drag from 
the fuselage and thrust from the pro- 
peller. The weight is another force and 
this always acts vertically downwards. 
Fig 4 shows all the main forces acting 
on a plane in flight. Notice that the 
plane is not actually moving exactly 
where it’s pointing and that the lift and 
drag forces are parallel and perpen- 
dicular to the direction of motion. 
Weight, of course, is always directed 
towards the ground and thrust is 
aligned in the direction of the propeller, 
i.e. the direction in which the plane is 
pointing. 

Adding all the forces together when 
a plane is in flight gives a resultant 
force, and this is not difficult to pro- 
gram if you have some basic knowl- 
edge of vector algebra. If the resultant 
is zero then the plane is in equilibrium, 
which doesn’t mean it’s at a standstill: 
it just means it will continue moving at 


Fig 3 Aeroplane Instruments 
Primary Instruments 


Altimeter: Height above sea level in feet 


the same speed and in the same direc- 
tion as previously. The plane can be in 
equilibrium when on the ground, ina 
climb or descent as well as in level 
flight. 

When a plane is flying, the amount 
of lift produced by the wings and tail 
depends on a number of factors. Firstly, 
the greater the speed of the plane and 
the larger the wings, the more lift is 
produced. The design of the wings 
plays a part — some wings produce 
more lift (but often more drag) and 
others produce less. Also, in general, 
the greater the angle the wings present 
to the airstream, the greater the lift — 
this angle is known as the angle of 
attack. Notice in Fig 4 that the tail has 
a negative angle of attack, resulting in 
negative ‘lift’ with the overall effect of 
pulling the tail downwards. Lift 
increases up to a certain angle of attack, 
at which the wing ‘stalls’, suddenly 
producing less lift, and usually with a 
fairly dramatic effect on the plane’s 
progress through the air. Most wings 
stall at around 15° to 18°. 

If there’s a resultant force from all 
the various component forces then the 
plane accelerates in that direction. This 
force can be used to plot the plane’s 
new position (the fourth step in the 
loop in Fig 2). Force and acceleration 
are related by the famous equation from 
Newtonian physics: P =.M x A, where 
P is the force on an object, M is its 
weight (or strictly speaking its mass) 
and A is the resulting acceleration. The 
acceleration changes the velocity of 
the plane, so if, for example, the resul- 
tant force is upwards and the plane is 
already climbing, it will start to climb 
more rapidly. If it is currently going 
down then it will start going down less 
quickly. The user will see the effects of 
the acceleration as both an increased 
value on the vertical speed indicator 


Vertical Speed Indicator: Rate of ascent or descent in feet per second 
Air Speed Indicator: Plane’s speed through the air in knots 
Directional Indicator: Plane’s heading in degrees magnetic 

Artificial Horizon: Plane’s pitch and roll attitude 

Turn Co-ordinator: Plane’s rate of turn in degrees per minute 


Engine Instruments 


Tachometer: Engine/propeller speed in revolutions per minute 


Oil Temperature Gauge: White/Green/Red 


Oil Pressure Gauge: White/Green/Red 
Ammeter, Suction 


Navigation Instruments 
Automatic Direction Finder (ADF) 
VHF Omni Range (VOR) 
Instrument Landing System (ILS) 
Global Positioning System (GPS) 
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and a faster-moving altimeter. 

The key features of a plane’s position 
are its position in space, expressed as 
latitude, longitude and height, along 
with its orientation (i.e. its angle of bank, 
pitch etc), velocity and weight. It is 
obvious that position, orientation and 
velocity must be taken into account in 
the simulation, but the weight of a plane 
greatly affects the way it flies, so for an 
accurate simulation it should also be 
considered. The weight changes as the 
flight progresses, as all planes burn fuel 
and thereby lose weight. Some planes 
can also lose weight by dropping bombs 
or passengers (normally attached to 
parachutes) and in certain weather con- 
ditions aeroplanes can also gain weight, 
by accumulating ice. 

Looking at the diagram in Fig 4 it can 
be seen that not all the forces on the 
plane act on the same point, and so in 
addition to the acceleration the forces 
can also twist the plane in any direction, 
which alters its orientation. In particu- 
lar, its changing orientation can affect 
the amount of lift generated by the wings 
and tail. If the combination of forces 
pitches the nose up then this will 
increase the angle of attack of the wings 
and create more lift, up to the stall angle 
at any rate. 

So much for the theory. Fig 5 shows 
the flight simulator in action. The ele- 
vators, ailerons and throttle are con- 
trolled by scroll bars, although to keep 
things simple, in this implementation 
we've kept the ailerons locked. This 
means that the plane cannot make turns. 
Fig 5 The Visual Basic Flight Simulator 
in action 


FS Minus One 
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Fig 4 The main forces acting on 
the aeroplane in flight 


Weight 


Tail lift 
Tail 
_ atime Angle of attack 


Fig 6 Program Code 

[Sub Form_Load () 

‘Initialise all constants 

‘SET SIMULATION ENVIRONMENT 

‘Make scroll bars fairly sensitive 
ailerons.Max = 20 

elevators.Max = 20 

throttle.Max = 20 

timerl.Enabled = False 

timerl.Interval = 1 ‘Thousandth of a sec 
‘Don’t bother clearing text boxes - they’11 soon get written in 


Wing 
Angle of attack 


‘SET ENVIRONMENT PLANE FLIES IN 

‘Definition of International Standard Atmosphere 
‘Set up air densities up to 11,000 metres 

atmosN = 12 

altitude(1) = 0!: density(1) = 1.225 

altitude(2) = 1000!: density(2) = 1.112 
altitude(3) = 2000!: density(3) = 1.007 


altitude(4) = 3000!: density(4) = .909 
altitude(5) = 4000!: density(5) = .819 
altitude(6) = 5000!: density(6) = .736 


altitude(7) = 6000!: density(7) = .66 

altitude(8) = 7000!: density(8) = .59 

altitude(9) = 8000!: density(9) = .526 

altitude(10) = 9000!: density(10) = .467 

altitude(11) 10000!: density(11) -414 

altitude(12) = 11000!: density(12) = .365 d 

‘After 11,000 metres we are in stratosphere ; 
magVar = 4 * DEGREES_TO RADIANS ‘4 degrees W in SE England 
‘Could have wind velocities at altitudes here too 
North(1) = 0!: North(2) = 1!: North(3) = 0! 

East(1) = 1!: East(2) = 0!: East(3) = 0! 

Up(1) = 0!: Up(2) = 0!: Up(3) = 1! 


‘SET UP THE PLANE 
‘Definition of a sort of Cessna FRA 150L here. 
‘Figures come from an Aerobat Flight Manual unless otherwise stated 


‘Distances in metres from firewall for (pitching) balance 
calculations 

emptyFirewall = 416! / 520! 

occsFirewall = 153! / 154! 

fuelFirewall = 65! / 61! 

bagsFirewall = 15! / 10! 

wingFirewall = .48! ‘Measured from drawing - May not be exact 
StabFirewall = 4! ‘Ditto 


‘Wing dimensions 

wingArea = 14.8! 

wingSpan = 10.11! 

wingChord = wingArea / wingSpan 
wingIncidence = 1! * DEGREES_TO RADIANS 
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‘Wing roll rates (a guess) 
wingVell = 50! * KNOTS _TO_ MPS: wingRolli = 20! * DEGREES_TO_RADIANS 
wingVel2 100! * KNOTS_TO MPS: wingRoll2 = 90! * DEGREES_TO RADIANS 


‘Characteristics for a Clark YH Aerofoil, without flaps 
‘(Aerobat actually has a NACA 2412 with Fowler flaps) 
wingN = 16 


wingAofA(1) = -4! * DEGREES_TO RADIANS: wingCL(1) = -.09: wingCD(1) = 
201: WingCP(1) = .8 

wingAofA(2) = -2! * DEGREES _TO RADIANS: wingCL(2) = .05: wingCD(2) = 
-009: WingCP(2) = .74 

wingAofA(3) = 0! * DEGREES_TO RADIANS: wingCL(3) = .2: wingCD(3) = 
-01: WingCP(3) = .4 


wingAofA(4) = 2! * DEGREES_TO_RADIANS: wingCL(4) 
-015: WingCP(4) = .32 
wingAofA(5) = 4! * DEGREES_TO_RADIANS: wingCL(5) 


-36: wingCD(4) 


-51: wingCD(5) 


-022: WingCP(5) = .295 

wingAofA(6) = 6! * DEGREES_TO_ RADIANS: wingCL(6) = .66: wingCD(6) = 
-033: WingCP(6) = .285 

wingAofA(7) = 8! * DEGREES_TO RADIANS: wingCL(7) = .8: wingCD(7) = 
-045: WingCR(7) =*.275 

wingAofA(8) = 10! * DEGREES_TO RADIANS: wingCL(8) = .94: wingCD(8) = 
-062: WingCP(8) = .27 


wingAofA(9) = 12! * DEGREES_TO RADIANS: wingCL(9) 
-083: WingCP(9) = .27 
wingAofA(10) = 14! * DEGREES_TO_RADIANS: wingCL(10) = 1.21: 


1.06: wingCD(9) = 


wingCD(10) = .103: WingCP(10) = .27 

wingAofA(11) = 16! * DEGREES_TO_RADIANS: wingCL(11) = 1.33: 
wingCD(11) = .125: WingCP(11) = .265 

wingAofA(12) = 18! * DEGREES_TO_RADIANS: wingCL(12) = 1.43: 
wingCD(12) = .146: WingCP(12) = .265 

wingAofA(13) = 19! * DEGREES_TO RADIANS: wingCL(13) = 1.36: 
wingCD(13) = .17: WingCP(13) = .275 

wingAofA(14) = 20! * DEGREES_TO_RADIANS: wingCL(14) = 1.26: 
wingCD(14) = .211: WingCP(14) = .29 

wingAofA(15) = 25! * DEGREES_TO_RADIANS: wingCL(15) = .97: wingCD(15) 
= .324: WingCP(15) = .33 


wingAofA(16) = 30! * DEGREES_TO_RADIANS: wingCL(16) = .81: wingCD(16) 
= .43: WingCP(16) = .37 


‘Stabiliser dimensions 

stabArea = 1.58 

stabSpan = 3.05 

stabChord = stabArea / stabSpan 


‘Unlike C150 no true elevators - instead an all-flying tail-plane 
(stabilator).. 

stabMinAng = -14 * DEGREES_TO_ RADIANS 

stabMaxAng = 8 * DEGREES_TO_RADIANS ‘So half way is -3 degrees as in 
c150 


‘Stabiliser Flight Characteristics. C150 has an unknown (to me) 
symmetrical section at -3 degrees incidence 

‘This is a NACA 009 symmetrical, Centre of Pressure values are a 
guess 


stabN = 15 

stabAofA(1) = -14! * DEGREES_TO RADIANS: stabCL(1) = -.87: stabCD(1) 
= .036: stabCP(1) = .25 

stabAofA(2) = -12! * DEGREES_TO RADIANS: stabCL(2) = -.89: stabCD(2) 
= .028: stabCP(2) = .25 

stabAofA(3) = -10! * DEGREES TO RADIANS: stabCL(3) = -.9: stabCD(3) = 
-O21:''stabCP (3): = 225 

stabAofA(4) = -8! * DEGREES _TO RADIANS: stabCL(4) = -.85: stabCD(4) = 
~018: stabcCP(4) = .25 

stabAofA(5) = -6! * DEGREES _ TO RADIANS: stabCL(5) = -.64: stabCD(5) = 
-014: stabCP(5) = .25 

stabAofA(6) = -4! * DEGREES_TO RADIANS: stabCL(6) = -.43: stabCD(6) = 
VOL1? (stabCr(6):-=..25 

stabAofA(7) = -2! * DEGREES_TO RADIANS: stabCL(7) = -.21: stabCD(7) = 


BLACK REFILL TWIN PACK £13.99 
BLACK REFILL TEN PACK £49.99 


COLOURED REFILL TWIN PACK £14.50 
(YELLOW, MAGENTA, CYAN, BURGUNDY, GOLD, GREEN) 


DESKJET CARTRIDGE (NEW) 
BJ10e CARTRIDGE (NEW) 
PRICES INCLUDE VAT AND P&P 


RECHARGED TONER 
CARTRIDGES 


HP/CANON SERIES I, II, Il 

EXTRA LIFE 

IBM 40lie HIGH YIELD 

HP Ill Si £60.00 + VAT 

HP/CANON SERIES I, I] (BLUE,RED) £60.00 + VAT 
INCLUDES DELIVERY 


£35.00 + VAT 
£70.00 + VAT 


ABSTRACT SOFTWARE 


Create your own games with... 
Games Programmers ‘Toolkit. 
3 * 3.5 HD Disks £8.00 


Shareware source code and lib’s for parallax 
scrolling. Soundblaster control, animation, 
sprites, joystick control and much more. 


Graphics “Programming. 

7 * 3.6 HD disks £15.00 

Source code and spec’s for GIF, PCX, BMP, 
TIFF and TAROS. Plus Paint packages and 


conversion routines. 
Or buy both for £20.00 
Guaranteed 14 day delivery. 
Catalogue available. 


Add £3.00 P&P. Cheques and P.O’s to: 
Abstract Software 
41 Christchurch Street 
Cambridge. CB1 1HT 
Tel - 0223 67692 
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The instruments are shown as text 
boxes and are updated as quickly as 
the processor can manage. The Restart 
button places the plane at a set position 
in a cruise at 2000 feet, and Exit has the 
usual effect. All the controls have fairly 
obvious names, and any special values 
that they may require are set up as part 
of the initialisation. 

The code behind the scenes is shown 
in Fig 6. Form_Load initialises all the 
constant values, describing the plane 
itself and the atmosphere that it flies in. 
Restart_Click puts it in a position where 
it’s cruising at 2000 feet, with controls 
set in a sensible fashion. Timer1_Timer 
is called as often as possible by a timer 
control in the form. It performs one 
iteration of the main loop, as shown in 
the flow diagram in Fig 2. The sub- 
routine calls in Timer1_Timer invoke 
the four main tasks in turn: Display, 
ControlInput, Calculate and Newposn. 

The simulator has a number of 
sophistications. It accurately models 
the flight characteristics of a well 
known wing section, with the correct 
dimensions, loading capacity and 
engine power for a Cessna Aerobat. It 
can compensate for decreasing air den- 
sity at altitude, and contains the cal- 
culations for weight and balance. 
However, it leaves a number of inter- 
esting exercises for the reader. Firstly, 
it’s prone to instability, especially when 
it has to cope with differing processor 
speeds. It’s easy for the simulated plane 
to pitch up (or down) beyond the stall 
angle, and then go out of control. 
Improved calculations of moments 
should fix this, or alternatively an auto- 
matic ‘stick-pusher’ which prevents 
stalling. The ability to use ailerons 
would be a big advantage as well, 
although this requires a more sophis- 
ticated use of vector algebra than the 
treatment we have given here. 

For any readers interested in 
exploring flight modelling further, the 
following books are highly recom- 
mended: 

Reims Cessna FRA150L Flight 
Manual (Reims Aviation) — technical 
specification of the Aerobat. 

Mechanics of Flight, AC Kermode 
(Longman Scientific and Technical) — 
a good elementary mathematical intro- 
duction to the dynamics of flight. 

Air Pilot’s Manual - The Aeroplane 
Technical, Trevor Thom (Airlife) — a 
layman’s description of how planes fly 
and the effects of controls. 

Vector Analysis, Murray Spiegel 
(Schaum Outline Series) — a mathe- 
matics book describing the use of vec- 
tors, making the implementation of a 
flight model a great deal easier. 


-01: stabCP(7) = .25 . 
stabAofA(8) = 0! * DEGREES_TO RADIANS: stabCL(8) = 0!: stabCD(8) = 
-009: stabCP(8) = .25 : 

stabAofA(9) = 2! * DEGREES_TO_RADIANS: stabCL(9) = .21: stabCD(9) = 
-01: stabCP(9) = .25 

stabAofA(10) = 4! * DEGREES_TO RADIANS: stabCL (10) 
= .011: stabCP(10) = .25 

stabAofA(11) = 6! * DEGREES TO RADIANS: stabCL (11) 
= .014: stabCP(11) = .25 

stabAofA(12) = 8! * DEGREES_TO_ RADIANS: stabCL(12) = .85: stabCp(12) 
= .018: stabCP(12) = .25 

stabAofA(13) = 10! * DEGREES_TO RADIANS: stabCL(13) = .9: stabCD(13) 
= .021: stabcP(13) +25 

stabAofA(14) = 12! * DEGREES_TO RADIANS: stabCL (14) 
= .028: stabCP(14) 225 

sStabAofA(15) = 14! * DEGREES _TO RADIANS: stabCL(15) = .87: stabCD(15) 
= .036: stabCP(15) = .25 

stabPitch = (1! / 2000!) * DEGREES_TO RADIANS ‘Pitch change per 
Newton metre per second 


-43: stabCD(10) 


-64: stabCD(11) 


-89: stabCD(14) 


‘Power comes from Rolls Royce 0-240-A with McCauley 1A135/BRM7150 
1.803 metre propeller 

minrpm = 800 ‘Idle speed on the ground 

Maxrpm = 2800 ‘Max permitted - not max possible! 

‘At 2,500 feet at 2,800 RPM engine consumes 33.7 litres per hour 
fuelrate = (33.7 * LITRES_TO KGS) / (2800! * 3600!) ‘Kgs per rpm pre 
second 

‘Engine is actually 130Hp = 97 kilowatts 

‘How this translates into actual delivered power from the prop I know 
not 

‘These figures give guessed max poss delivered power (watts) at sea 
level 

engAlt1 = 0: engN = 7 

engveli(1) = 50 * KNOTS_TO MPS: engPow1(1) = 41000 


engVel1(2) = 60 * KNOTS_TO MPS: engPow1(2) = 57000 
engVel1(3) = 70 * KNOTS_TO_ MPS: engPow1(3) = 65000 
engVeli(4) = 80 * KNOTS_TO_ MPS: engPow1(4) = 70000 
engvel1(5) = 90 * KNOTS_TO_ MPS: engPow1(5) = 75000 


engVel1(6) = 100 * KNOTS_TO_MPS: engPow1(6) = 80000 

engVel1(7) = 110 * KNOTS_TO_MPS: engPow1(7) = 82000 

‘These figures give guessed max poss delivered power (watts) at 
service ceiling 14,000 feet 


engAlt2 = 14000 * FEET _TO METRES . 

engvel2(1) = 50 * KNOTS_TO MPS: engPow1(1) = 20000 
engVel2(2) = 60 * KNOTS_TO MPS: engPow1(2) = 27000 
engVel2(3) = 70 * KNOTS_TO MPS: engPow1(3) = 33000 
engVel2(4) = 80 * KNOTS_TO MPS: engPow1(4) = 35000 
engvel2(5) = 90 * KNOTS_TO_ MPS: engPow1(5) = 37000 


engVel2 (6) 


100 * KNOTS_TO_MPS: engPow1(6) = 40000 
engVel2(7) = 110 * KNOTS_TO_MPS: engPow1(7) = 41000 
CFD = .1’A guess! 

Call restart_Click 

timerl.Enabled = True 

End Sub 


Sub Calculate () 

Dim rho As Single, speed As Single, q As Single 

Dim theta As Single, ang As Single, thi As Single 

Dim Cl As Single, CD As Single, CP As Single 

Dim wingLift As Single, wingDrag As Single, wingCPFirewall As Single 
Dim stabLift As Single, stabDrag As Single, stabCPFirewall As Single 
Dim fuseDrag As Single 

Dim totThrust As Single, totDrag As Single, totWeight As Single 

Dim totLift As Single 

‘Calculate all forces and moments acting on plane 

rho = AirDensity(P(3) ) 

speed = VMag(V()) : 

q= .5 * rho * speed * speed ‘used in several equations below 

‘Calc Angle of attack relative to thrust line T using 
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ey eh IT| IVI cos theta (and |TI = 1) 

theta = arccos(VDot(T(), V()) / VMag(V())) 

“Calc Wing Coeffs 

ang = theta + wingIncidence 

Call Coeffs(ang, wingN, wingAofA(), wingCL(), wingCD(), WingCP(), Cl, 
CD, CP) 

‘Calc Lift, Drag and Moment on Wing 

wingLift = q * Cl * wingArea 

wingDrag = q * CD * wingArea 

wingCPFirewall = (wingFirewall + wingChord * (1! - CP)) 

‘Calc Elev Coeffs 

ang = theta + elevangle 

Call Coeffs(ang, stabN, stabAofA(), stabCL(), stabCD(), stabCP(), Cl, 
CD, CP) 

‘Calc Lift, Drag and Moment on Stabiliser 

stabLift = q * Cl * stabArea 

stabDrag = q * CD * stabArea 

stabCPFirewall = (stabFirewall + stabChord * (1! - CP)) 

‘Calc Drag on Fuselage 

fuseDrag = q * CFD 

‘Total Litt 

totLift = wingLift + stabLift 

‘Total Drag 

totDrag = wingDrag + stabDrag + fuseDrag 

‘Total thrust 

totThrust = Thrust(P(3), engRpm, speed) 

‘and weight 

totWeight = (empty + occs + fuelwt + bags) * G 

‘Angle between Velocity V()and Up() 

‘given by V. Up = IVI |[Up| cos thi (and |Up| = 1) 

thi = arccos(VDot(V(), Up()) / VMag(v())) 

‘Now the resultant forces 

Res(1) = 0! ‘Doesn’t move sideways 

‘Horizontal 

Res(2) = totThrust * Sin(thi - theta) - totLift * Cos(thi) - totDrag 
* Sin(thi) 

‘Vertical 

Res(3) = totThrust * Cos(thi - theta) + totLift * Sin(thi) - totDrag 
* Cos(thi) - totWeight 

‘Work out moments using baseline in dirn of V assume thrust and 
fuselage drag act thru centre of firewall 

totMom = (wingLift * wingCPFirewall + stabLift * stabCPFirewall) * 
Cos (theta) 

totMom = totMom + (wingDrag * wingCPFirewall + stabDrag * 
stabCPFirewall) * Sin(theta) 

totMom = totMom - (empty * G * emptyFirewall - occs * G * 
occsFirewall - fuelwt * G * fuelFirewall - bags * G * bagsFirewall) * 
Sin(thi) 

totMom = totMom + (empty * G * emptyFirewall - occs * G * 
occsFirewall - fuelwt * G * fuelFirewall - bags * G * bagsFirewall) * 
Cos (thi) 

‘pitchchange = stabPitch * totMom * duration 

‘Lots of stability problems with the above, so.... 

PitchChange = (theta + elevangle) 

‘End Function 


Sub Coeffs (AofA As Single, N As Integer, Angs() As Single, CLift() 
As Single, CDrag() As Single, CPress() As Single, Cl As Single, CD As 
Single, CP As Single) 

‘Returns Coefficient of lift, drag and Centre of Pressure 

‘for a given Angle of Attack (in radians) 

Dim i As Integer 


Lee 

Do While i < N And Angs(i) < AofA 
feo 1 

Loop 


Cl = LinInterp(AofA, Angs(i - 1), CLift(i - 1), Angs(i), CLift(i)) 
CD = LinInterp(AofA, Angs(i - 1), CDrag(i - 1), Angs(i), CDrag(i)) 
CP = LinInterp(AofA, Angs(i - 1), CPress(i - 1), Angs(i), CPress(i)) 
End Sub 


Barcode Solutions. 
Any PC, Any application, 
Windows or DOS. 


‘Most applications can benefit from faster, error free data entry by 
barcode. Altek systems need no complicated set up, no thick manuals to 
read, simply plug and go. Our package deals are complete and include 
everything to get started including software to print barcodes on HP. 
LaserJet and Epson Dot Matrix printers, All types of barcode scanner are 
available. 


We were first to introduce Barcodes to PCs. Many public sector 
and private companies both large and small now benefit from 
‘Altek systems. Larger systems for factory or offices are 
javailable. 


|Altek systems are effective, simple to operate, reliable 
and reasonably priced. 


Ask for further details 
for contact us with 
jan outline of 
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Enterprise House 

44-46 Terrace Road 
Walton-on-Thames, Surrey 
KT12 2SD 

Tel: (0932) 244110 

Fax No. (0932) 227279 


The Institution is Britain's leading specialised professional 
body for systems analysts and computer programmers. It 
now has some 3,000 members in 60 different countries. 


The purpose of the Institution is to assist members to 
advance their careers, and to secure public recognition for 
the valuable work they do. Members are required to 
maintain high professional standards. 


Apart from the prestige which Institution membership 
bestows, benefits include newsletters and other specialised 
publications, plus access to specialists who can assist with a 
wide range of educational, technical and business problems. 


The Institution of Analysts and Programmers is a non profit- 

making company owned by its members and governed by an 

elected Council. Applications for membership are welcome. 
* 


Tel: 081-567 2118 Fax: 081-567 4379 


* 
THE INSTITUTION OF 
ANALYSTS & PROGRAMMERS 
Charles House, 36 Culmington Road 
London W13 9NH, England 
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HANDS ON - COMPUTER ANSWERS \ 


Do you have a computing 
problem, or think you could 
help other readers out? 
Contact our resident expert, 
Frank Leonhardt. 


Procomm problem 


Some time ago I fitted an internal 
fax/comms modem in my 386 PC to 
replace a comms-only modem. The new 
modem was packaged with Winfax 
Lite. I have used Procomm ever since 
I started using modems over a year ago, 
and still use it with the new modem. 
The problem I have is that if the 
modem is used in fax mode and in the 
same session thereafter in comms 
mode, the modem still thinks it’s a fax 
and consequently tries to hand-shake 
with the modem at the other end in fax 
mode. There is no problem the other 
way around. The only way to overcome 
the problem is to re-boot the PC with 
the three-finger salute (Ctrl+Alt+Del). 
I suspect that the problem lies with 
Procomm not resetting the modem cor- 
rectly to the comms mode. Is there a 
way to overcome this (i.e. is there a 
routine which can be added to 
Procomm)? 
Stan Smith, Preston 


If it’s any consolation, this is a common 
problem. Manufacturers of fax software 
tend to assume that you will only be 
using your modem for fax in the future 
and therefore don’t reset it to data mode 
when they have finished. 

The command which should set 
fax-modems back to data is 
AT+FCLASS=0 (check this is correct 
before printing!), and you could try 
adding it to Procomm’s modem ini- 
tialisation string. I haven’t always had 
complete success with it. The modem 
supplier ought to tell you how to reset 
the modem without resetting the whole 
machine, but I expect you’ve already 
tried this approach and met the same 
indifference and confusion as I did. 


Gold crash 


I’m having repeated crashes with Gold 
Disk Animation Works. My machine 
is a Ti’Ko 486DX running at 33MHz 
with 8Mb RAM. The video card is an 
Orchid Fahrenheit 1280. 

From the reams of log generated by 
good old Dr Watson, it would appear 
that the program is clashing with the 
video driver. The program is only sta- 
ble if 1 run Windows at standard VGA 
resolution. 


The people at Gold Disk have tried 
their best to sort it out, but with little 
success. Their advice is to work at VGA 
resolution! Not very appealing — too 
few colours. 

I wrote to Orchid and sent them a 
few pages of log so that they might be 
able to solve the glitch. Their response 
was to send me the latest video driver 
5.11. Did it work? I’m afraid not. The 
program keeps on crashing, especially 
when I try to use the cell editor. 

It seems nobody is prepared to accept 
responsibility. Gold Disk say it’s Orchid’s 
problem and Orchid say it’s up to Gold 
Disk. Meanwhile I’m in the middle, 
with a very flaky piece of software. 

Diane Colover, W Yorks 


I wish I understood just what it is about 
Windows’ video drivers which gives 
hardware manufacturers so much trou- 
ble. Your problem is typical, in that 
you have a piece of application soft- 
ware which works fine using the stan- 
dard VGA.DRV and which crashes 
using someone else’s. The thing to 
remember about VGA.DRV is that it 
has been around a long time and has 
therefore been tested with an awful lot 
of software. The same isn’t true of any 
other drivers. 

Applications which make use of DOS 
extenders seem to have the most prob- 
lems. Even Microsoft’s own visual C++ 
1.5 comes to grief with a standard IBM 
XGA-2 display, though it seems to work 
with most others. 

If you can’t stand going back to 
standard VGA, it’s probably worth 
trying the Méicrosoft-supplied 
SGVA256.DRV or SUPERVGA.DRV 
instead. If they don’t work then you'll 


Particularly since MSDOS 6.0 was 
launched, with DoubleSpace disk 
compression built in, some 
manufacturers and dealers have started 
listing their hard disks, tape streamers 
and so forth in terms of the amount of 
‘storage they offer. When you read this 
small print, you often find a statement 
along the lines of ‘supplied with free 
software making the drive capacity 
XMb’, where X is the genuine capacity 
multiplied by two. | 

Using a software compression driver 
like Stacker, SuperStor, DoubleSpace. 
or the current offering in DRDOS does 


_ watch out. 


al PERSONAL COMPUTERWORLD SEPTEMBER 1994 


just have to wait for Gold Disk and 
Orchid to get their acts together. 

It’s very difficult to say who is to 
blame in a situation like this. These 
types of problem are usually down to 
the application programmers taking a 
liberty with the display driver which 
just happens to work with VGA.SYS. 
Whoever wrote the driver for the Orchid 
would have made sure that it accepted 
the same function calls as VGA.DRV, 
and done the same things regarding 
memory management. They will not 
necessarily have copied any function- 
ality which was not explicitly detailed 
in Microsoft’s documentation. 

So you could say that it’s Orchid’s 
fault for not copying VGA.DRV’s fea- 
tures exactly, or you could say it’s Gold 
Disk’s for using something in VGA.DRV 
which wasn’t properly documented. 
Take your pick. 

It would save a lot of time and trou- 
ble if driver and application program- 
mers got together in order to find out 
exactly what's going on. I suspect that 
all these incompatibility problems stem 
from a few simple snags which it must 
be possible to find, document and avoid 
in the future. Meanwhile, if you suffer 
from mysterious software crashes under 
Windows, try going back to VGA.DRV. 


Failed floppy 


Iam having problems with the floppy 
disk drive in my computer, a 33MHz 
486 PC with a Cyrix chip. Copying large 
files to it gives the error ‘Data error 
writing to drive A:’. Choosing ‘Fail’ at 
the ‘Abort, Retry, Fail?’ prompt gives 
‘Fail on INT 24’. 

Changing the disk drive, controller 


Buyers Beware: Double-Sized Disk Drives 


compression softwar 
start, you may not t 
compatible PC. It m 1 
unacceptable reliability, perform 
and memory penalties, or y 

may be unsuitable 

other reasons. _ 


